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Abstract 

We give the connections among the Fekete sets, the zeros of orthogonal 
polynomials, l(TO)-normal point systems, and the nodes of a stable and 
most economical interpolatory process via the Fejer contants. Finally the 
convergence of a weighted Griinwald interpolation is proved. 



1 Introduction 

L. Fejer introduced the so-called Hermite-Fejer interpolatory process, and in 
1934 he gave the definition of normal- and (p-normal system of nodes for which 
the Hermite-Fejer interpolation is a positive interpolatory process. The sur- 
prising nice convergence properties of Lagrange, Hermite and Hermite-Fejer 
operators on g-normal systems were proved by L. Fejer, G. Griinwald, etc. On 
the other hand the experiences in electrostatics ensure a system of nodes: the 
Fekete set, which has uniform distribution in some sense, so it must be a good 
set for interpolation. The system of zeros of orthogonal polynomials has very 
similar properties, as it it well-known. From another point of view, Egervary 
and Turan asked, that is it possible to find an interpolatory process, and a sys- 
tem of nodes together, such that the interpolatory polynomial has the minimal 
degree, and the operator has the minimal norm. The above-mentioned point 
systems can be a suitable system of nodes for an interpolatory process in general 
sense and also with respect to the Egervary- Turan problem. 

The primary aim of this note to revisit the connections among that sets of 
nodes, and interpolatory problems investigated e.g. in [2], [3], [4], [5], [6], [10]. 
In the next section, we summarize and reformulate these results, and complete 
them, when the original statement proved only in classical cases. It will be 
pointed out, that in these equivalences the so-called Fejer constants (see(3)) 
play the key role, that is the characterization of this special system of nodes is 
ensured by the Fejer constants. 

As an application of the results of the second section, in the third section 
we prove a convergence theorem on Griinwald interpolatory process on the real 
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line for Freud-type weights. As it turned out, giving the weighted Fekete sets 
with respect to a fixed weight is difficult. (However, there are several methods 
of giving approximating Fekete sets.) The zeros of orthogonal polynomials are 
Fekete sets for some varying weights. Unfortunately these varying weights tend 
to zero locally uniformly, so interpolation on Fekete sets in this sense gives only 
trivial (convergent) processes. The investigation of these weights at infinity 
leads to define a weighted Griinwald operator (see (11)), which has rather nice 
convergence properties. Comparing this result with the previous ones of [5], 
[13], it turns out that the convergence is valid here for a wider function class. 



2 Connections 

At first we give the definition of classes of weights in question. 

Definition 1 Let Y, d <C is a closed set. w is quasi-admissible on E, if it is non- 
negative, upper semi- continuous, and if T. is unbounded, lim|z|^oo |z|w(z) — 0. 

It is admissible, if cap {z g S : w{z) > 0} > 0. Let us call an admissible weight 
as approximating on (a, b) C M, if it has finite moments, it is twice differ entiable 
and (log (:^)) " > on (a, b), and if a is finite, then \mvx^a+ -^i^ = 0, and if 
b is finite, then Imix^b- — 0- 

Definition 2 \1 If llL 1 Let w be a quasi- admissible weight on a closed set S C 
C. Then 3^„ are called n-th weighted Fekete sets associated with w, if the supre- 
mum below is attained at the set Jn.w = {^i, ■ ■ ■ , Xn}- 



dn,w — sup drfi^^{^Z\, . . . , Zji'j 

zi,...,z„eii 



= sup TT \zi - Zj\w{z^)w(^Z.j) \ (1) 

Usually these points arc not unique, but in one dimension by some restric- 
tions on the weight, uniqeness can be proved. In the classical, unweighted case 
on [—1, 1], the result is proved by Popoviciu (cf. [U] Ch 6.7 p. 139., and the 
reference therein). In weighted case, after some restrictions on the weight a 
representation of Fekete points was given by M. E. H. Ismail ([F Thms. 2.1, 
2.4), wich ensures the unicity of the Fekete sets as well. In the foUowings the 
one-dimensional case will be investigated. 

Now let us deal with the weighted Lagrange interpolatory polynomials on 
a system of nodes X = {xkn,k = l,...,n:n e Nj. Let lk(x) ~ ,, ^[f^ r, 

where u;(x) = nfc=i('^~ -^fe) (denoting Xk — Xk^m k — 1, . . . ,n) the fundamental 
polynomials of the Lagrange interpolation, and let w{x) = e"'^^^-' be an ap- 
proximating weight. The properties of L^ ^ x{x) — L^.tuix) — w(x) ^ \ will 
be investigated. It is clear, that Lk.w{xk) — 1, that is the sup-norm of this 
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weighted polynomial is at least 1. If this sup-norm is equal to one, then Lfe^^ 
has a maximum at the point Xfc, that is 



= wyx) 

which ensures that 



w{xk) \ w{xk) lk{xk) 



w{x)^^[-Q{xk) + ^{xk)\=Q, (2) 



Ck,w ■^Ck,w,x ^'^{xk) + ^{xk) =0- (3) 



This is the case, when X is a Fekete set with respect to w^c^-i) ^ namely 
ni<i<.. (^{x - xi)'^w'^(^{x)w"^^^{xi)j 



L 



k,w,X 



l\i<i<^ ({Xk - 2:/)2w2(„-i) (a;fc)i(;2(„-i) (j;,)") 

l^k \ / 

ni<«<.<" ({xi - Xj)^w^i^ {x,)w^i^ {Xj)) 
n i<^<3<r. Uxi - x-j Y-w"^^ (a;i)w^(^ (xj)) 



< 1, 



because in the denominator appears d ^ \ . 

It will turn out in the followings, that the behavior of the constants Ck,w as 
an indicator, shows the properties of the point systems, interpolatory systems 
and operators. Emphasizing the importence of these constants, let us call them 
as "Fejer constants". 

Following carefully the proof of the above mentioned theorem of Ismail 
([S],Thm. 2.1), we get the following 

Proposition 1 Let w be an approximating weight on an interval (a, 6). Then 
rf"*'" ^\ (zi, . . . , Zn) attains its maximum on (a, b) at a unique set 3^ i , 

for which the following characterization is valid. 

J _L = {xi, . . . ,x„} if and only if C^^ = 0, fc = l,...,n. (4) 

At first we have to note here, that finite moments are not necessary in this 
statement. According to Ismail [S], the proof of this theorem is the following: 
taking the partial derivatives of log d"^" ^\ , it turns out, that 

^logd"*-" "^! (xi,...,x„) = j = if and only if Ck.w =0, k = 

1, . . . , n. Computing the Hessian, it can be seen, that — iJ is always positive 
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definite, so reealling the boundary condition on w, we get that the maximum- 
set is unique, that is it is the unique solution of the equation system: Ck,w = 
0, A: = 1, . . . , n. Independently of the previous chain of ideas, an elementary 
proof on unicity can be given. 

Proposition 2 Let w be an admissible, continuous weight on R such that log ^ 
is convex. Then the associated weighted Fekete sets are unique. 

Proof: Contrary, let {xi}f^i and {yi}f^i are Fekete points with respect to w 
enumerated in increasing order, and let Zi = . Then because of the ordering 
of the points, and the log-convexity of the weight, by the arithmetic-geometric 
mean inequality 



\zi - Zj\w{Zi)w{zj) 



ix,-Xj) + ivi-V,) 



yt\ / Xj + yj 



\xi + \yi -yj\ f x,+ y, \ (xj+ y^ 
^— —w ] w ' 



> ^\x.^ - Xj\^ \y^ - y^y w{x^)^/ w{yi)^ w{xj)^ w{yj), 

where the inequality is an equality if and only if Xi = yi for all indices, wich 
establishes the uniqueness. 

For special weights, the Fekete sets are the zeros of some orthogonal polyno- 
mials (cf. [5], Before setting the precise statement we need some definitions. 

Definition 3 Let w = e^'^ be an approximating weight on {a,b). Let 

An{x) = gn f plJt)wit)9M-^dt, (5) 

J a IX 

where Pn,w — JnX^^ + . ■ . is the n*'' orthonormal polynomial with respect to w, 
and Qn = 

Now we can define our weights: 

Definition 4 Let w be as in the previous definition. 

In the following investigations the constant g„ has not any role, but it will 
come into the picture inconnection with a convergence theorem in the next 
section. Let us see some examples on ^l^^^l ([5]), which in classical cases are 
different only in normalization from the weights wi (jlOj). for which the deriva- 
tives of Pn.w-s are orthogonal : 
Example: 



4 



(1) If w = e-^', 

Qn 

that is Wn — \w independently of n and x, and here Wi — w — 2w„. 

(2) If w = a;"e-^, 

An{x) ^ 1 

Qn x' 

that is Wn = x^^^e^^ = xw independently of n, and here wi — Wn- 

(3) If w= (l-a;)"(l + a;)'^, 

An{x) _ a + l3 + l + 2n 
Qn 1 - a;2 

that is Wn = a+f}+i+2n (1 ~ ^)"^^ (1 + xY^^ , and here Wi — {a + (i + l + 2n)wn- 

(4) If w = e-^', 

■An{x) r,/ 2 I 2 I 2 \ 

— = 2(a; + Qn + Qn+l). 

Qn 

that is Wn = 2(^2+^J+^«^- 

From another point of view Wn has also an importance. Denoting by Pn^/Wn = 
Zm it satisfies the following differential equation with some $„ (cf. [5], Th. 3.6.): 

Zn{x)+^n{x)Zn{x) =0 (7) 



In the next statement we reformulate the results of Ismail, Rutka and Smarzewski 

(cf. [5], m)- 

Proposition 3 Let Wn be as in the definitions above, and let us assume that 
Wn is an approximating weight. Then 

Ck,w„ — 0, k = 1, . . . ,n if and only if {xk} the zeros of pn.iv (8) 

The proof of this statement depends on the differential equation of orthog- 
onal polynomials. The equation system on Ck-s means that the differential 
equation fulfils at the points Xk, k — I, . . . ,n. In the classical cases, it is a 
Sturm-Liouville equation, that is there are polynomials of degree n in the dif- 
ferential equation, which is realized at n points. In general cases unicity is 
used. 

Normal and g-normal point systems were introduced on [—1, 1] by L. Fejer 
in 1934 ( 2 ). The weighted analogon of this definition was given in [3j. The 
original aim of these definitions was assuring the positivity of the Hermite-Fejer 
interpolatory operator. The limit case, when q — I was investigated on the 
weighted real line in [3]. Here this last definition is cited only. 
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Definition 5 Let w be an approximating weight on (a, &). A system of nodes 
X — {xk,m k — I, . . . ,n;n G N} is l^w) -normal, if there is an L > 1 such that 

I < La„, (9) 

where a„ is the M-R-S number, and 

where lk{x)-s are the fundamental polynomials of the Lagrange interpolation. 

In this definition the kernel function of the Griinwald operator appears. 
Here we will follow the notations of |5] and [13] , that is the weighted Griinwald 
operator on the nodes {xk\^^i with respect to an / is 

n 

w{x)Yr,{f,x)=w{x)Y,ll(x)f{xk) (11) 
k=l 

Mostly the boundedness of the operator-norm ensures the convergence of the 
interpolatory process. The boundedness by one, is a very special criterium. 
This is the case for instance, when the reciprocal of the weight function has 
non-negative even derivatives, and the Griinwald operator coincides with the 
Hermite-Fejer one. Also on this chain of ideas the Fejer constants play the 
key role. More precisely, with the notations above, the weighted Hermite in- 
terpolatory polinomial (with some weight w) of a differentiable function can be 
expressed as (cf. @]) 

and the corresponding weighted Hermite-Fejer operator is 

u;(x)g„,^(/,x) =u;(^)f] ^^'^^•'"'^;~'^^^^^'^''\ ,M(:.fc), (13) 

which coincides with the weighted function at the nodes {xk}^^i^ and which 
has zero derivatives at the nodes. Furthermore by the definition of the Fejer 
constants, Hn^w{f ,x) is the (unweighted) Hermite interpolatory polynomial of 
— . So when the Feier constants are zero 

w J 

Yn,w{x) ■■= w{x)Yn [ —,x] = w{x) — = w{x)Hn,w{ — , x) 
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^ ^ H-c^A^^ xkMi^) ^ i, f 1) , . (14) 

is the Hermite interpolatory polynomial of ^ with respect to the nodes: {xk}2=i- 
So the following connections are established. 

Proposition 4 Let w be a weight as above. 

If a system of nodes {x^} is l[w)-normal, then Ck.w — fc = 1, . . . , n. 

On the other hand, let us suppose further, that (;jj)^" > on \x\ < La„. 
Now 

if Ck,w = 0, fc = 1, . . . , n then the system of nodes is l{w) -normal. 

Proof: If {xfc} is l(z«)-normal, then z/;(x) < 1, fc = 1, . . . , n (see (10)), so 
Ck,w =0, fc = 1, . . . , n. According to (14), by the error formula of the Hermite 
interpolation, it is clear, that 1 — w(x) ^^Lj^ u]{xl) — ^' '^hen (:^)^^"'' > on 

\x\ < LOn. 

The Egervary-Turan interpolatory problem (cf. |10| . and the references 
therein) is to find an interpolatory process of lowest degree, and of smallest 
norm. Below we denote by ^^(2;) each polynomial of arbitrary degree for which 
lk{xi) = 6ki, i = 1, . . . , n. 

Definition 6 Let w be as in Definition 5. The interpolatory system of polyno- 
mials lk{x), k — 1, . . . ,n is w-stable on {a,b) if for all yi, . . . ,yn > 

" ik (x) 

0<w{x)y^ — — r2/fc<maxyfc, x e {a,b). (15) 

^ W[Xk) k 

A w-stable interpolatory system on (a, b) is most economical, if 



E 

fc=i 



deg {ik{x)) (16) 



is minimal. 



Let us remark that if the weight function tends to zero quickly at the bound- 
ary points of the fundamental interval, then the w-stability of the Griinwald 
operator coincides with the l(w)-normality of the nodes. It is proved for all the 
classical weights (cf. [10], Thm. 2.3), that an interpolatory system is Wn-stable 
and most economical, if and only if it is the Griinwald operator on the zeros of 
Pn,w. From the previous investigations, similarly to the classical cases, we can 
state the parallel theorem for general weights. 

Let us denote by 
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Proposition 5 Let w be an approximating weight on an interval (a, &). 
If In{x) is w-stable and most economical, then 

Ck,^=0, fc = l,...,n. (17) 

Let us assume further that (^)^^"'' > on \x\ < Lon- 
If Ck,w = 0, fc = 1, . . . , n, then 

In{x) = Yn^w{x) is w-stablc and most economical (18) 

Proof: As it was pointed out eg. in [10 , if an interpolatory process /„(x) is 
w-stable and most economical, it must be the Griinwald operator, because by 
the positivity of the operator, lk{x) has zeros at the points Xi,i — 1, . . . ,n,i k 

of even multiplicity, that is X]fc=i 'i^S {jkix)j > 2n{n — 1). It is realized by 
Yn^w As it was shown in Statement 4, if Yn,w has maxima at x^-s then C^-s 
are zero. The opposite direction is also follows from Statement 4. 

Finally enumerating the properties discussed above, we can summarize these 
results as it follows. 



(A) Ck,n,^0, fc = l,...,n (A') Ck,n,„^0, fc = l,...,n 

(B) J _i = {xi,...,Xn} (B') J _i ={xi,...,Xn} 

(C) p„,w{xk)^0, k = l,...,n 

(D) Yn_w{x) is w-stable and most economical 
(D ) Yn_w„ (x) is Wn-stable and most economical 

(E) {xi,...,Xn} is l(w)-normal (E ) {xi,...,Xn} is l(u;„)-normal 

Through the equivalence of all the above mentioned properties with property 

(A) (or (A )), that is the Fejer constants are zero, one can get 

Corollary: 

Let w be an admissible, approximating weight on an interval (a, 6). If 

\ / \ (2n) 

ii;) > on (a, 6), then (A),(B),(D),(E) are equivalent, and if f > 

on (a, 6), then (A'), (B'), (C), (D'), (E') are equivalent. 

We have to show an example on the second assumption . 
Example: Let 

m 

Q(x)=^dfcx2^ dk>0, fc = l,...,m, (19) 

and let w{x) — e~'^^'^'. For these special Freud-type weig hts M- > 
on R for all n E N . According to the Leibniz rule it is enough to show that 
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it) ii^f"~'^ > for j = 1, . . . , 2n. Because 



2fe-2 



where 6;-s are positive, taking into consideration that w is an even weight func- 
tion, (and so Pn{w) is also even), one can see that 



III, ^rv — „ 



is a polynomial of x with nonnegative coefficients, and all the exponents of this 
polynomial arc; civen if j is even and are odd if j is odd. By a simple induction 
one can see that 

1 \ (^'^ 

-- =p(j,x)eQ(-\ 
^w{x) J 

where p{j, x) is a polynomial having the same properties as the previous one. 
Because j and 2n — j have the same parity. 



\ (2n) 



— ) {x) = p{x)e'^^''\ 

where p{x) is a polynomial with even exponents and positive coefficients, so it 
is positive on the real line for all n G N. 

Finally we have to remark that the assumption (^)*'^"'' > seems to be 
assymetric, and it is necessary only because of the method of the proof by 
Hermite interpolation. The question that can it be weakened or not, is unsolved 
yet. 



3 Interpolation 

In this section, let w = be a three times continuously differentiable Freud 
weight on M, that is we suppose that Q is even, Q > on (0, 00), and for some 

A,B >2; A< ^^Q.(^)^ < B on (0, 00), moreover there is a constant c such that 
for every \x\ > 1, ^^^^-^ < c. By these assumptions there is a d > 1 such 
that Q"{x) > ^-{B-i)^ -c{B-i) ^ ^^^^ > ^ 

we can define 

Definition 7 With d> 1 given above, let 



w 



w{x), \x\ < 1 

i^) = { ^^Wq^'N^'^ (20) 
twice continuously differentiable, elsewhere 
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Furthermore we assume that log i has positive and continuous first and 
second derivatives on (0, cxd). 

Let us remark at first that Q' {1} < Q' (d) + (^^^^ (d), because Q' 

is increasing, and the second member of the right-hand side is positive, when 
A > 2. That is a suitable connection can be defined between the two parts of 
log 4. 

As usually we define 

Definition 8 

C^^{feC{R)\ lim {fw){x)^0 

Let Yn{f, •) be as in (10), the Griinwald operator on the zeros of Pn,w Now 
we have the following 

Theorem 1 Let f e Then 

lim ||(y„(/)-/)zD|j =0 (21) 

n— ^-oo 

Comparing this theorem with Cor. 2. of 13;, we can see, that we have two 
different weights in this theorem, but when A > 2, then the function class is 
wider here, that is the fuctions can grow more quickly at infinity. 

The previous definition of the weight was inspirated by the next lemma. In- 
vestigating the weights Wn from the previous section, it turns out, that however 
Wn tends to zero locally uniformly when n tends to infinity, the behavior of Wn-s 
are the same at infinity. It means, that the Griinwald operator on Fekete points 
with respect to the varying weights w„ has trivial convergence properties, but 
it allows to find a non-trivial process, as it is given in the theorem. 

The following estimation of An is valid. 



Lemma 1 Let w be as above, and let A > 2. Let Lq be a constant such that 
■^a„ > a2n+[A~i]+i- For every L > Lq 

Qr, ^[^, zf \x\ >Lan ' 
where the constants in " ^ " depend only on L, but they are independent of n. 
Proof: 

At first we have to note that such an Lq exists by [7] 5.9. The first line 
of the inequality is proved by H. N. Mhaskar ([9], Prop. 3.7.). To prove the 
second line we have to divide the integral to some parts. Since An is even we 
can choose x > ia„. 

^ ^ - ' plit)w{t)-^^ — ^-^dt+ / {.)dt^Li+L2 

|t|<f J\t\>§ 
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If |i| ^ f J using that Q is convex on R, and estimating the denominator as 
%<\t-x\< |a;, 

2 (2^-^-1) Q'jx) ^ 2 Q'{x)-Q' {%) ^ Q'{t)-Q\x) 
3 2-4-1 a; - 3 x " t-x 

^ ^Q\t)\+Q'{x) ^ ^Q'{x)^ 
where in the first inequahty we used the properties of Q , of [7] 5.3. So 

h / P„w ^ / PnW , (23) 

X J\t\<§ 3; a; 

where in the second "tilde", the lower estimation fulfils by '9^,(2.6), say. 
For 1^1 > ^a„ > 4a„ 

l2= [ Pl{t)wit)^^^\^^dt= (.)+ (•)=/3+/4 (24) 



t — X J l'l>5 



-ti <i 



By the properties of Freud weights, and by [9 ,(2.6), 

h = . pl{t)w{t)Q" {^{x,t))dt 



- ''^i^ / Pnii')'^ii)dt < cie^^^"^^ (25) 
■ , ,Q'{\t\) + Q'(x) , , ,Q'(t) - g'(a;) , 

p^(t)7z;(t) ^ il^ / P^(^)u'(t) ^ ^ ; U t 

'\l{t)w{t)9Ji}-RMdt+ r{.)^h + I,+I, + I, (26) 



x + l ^ X J2x 



According to [T, 5.2, 



Q {x) 



' X 



According to 2.7, 



pl{t)w{t)t^^-'^^+^dt < cie-^^" / pl{t)wit)t^^-^^+^dt 

J lt|<a2„+[yi_i] + i 
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Because A> 2, ^-j^ is increasing, so by [T 5.2 and 5.9 

pl{t)w{t)t^^-^^+^dt < / pl{t)w{t)t^dt 

t|<02n+[A-l] + l ^ •^I*l<a2n + [A-1] + 1 



that is 



h < C3e-=-"^ (27) 



Ift€ — l) we can write that x =^ Xt, where 1 < A < 2, that is according 
to [7], 5.3, recalhng that B > A > 2, ^'^"^if < ^^^1^' < c{B)^ < 
c{B) '^ ^\ where in the last step we used that is increasing. So 

I^<c^M. r \l{t)w{t)dt<cie-^-^^^ (28) 

X J ^ X 

2 

As in the previous case, when x + 1 <t <2x, ^ < ^^^^^^^ so 

/7 < c{B)^^ r pl{t)w{t)dt < cie-^^"^^ (29) 

X Jx+l ^ 

Because t — x > | in /§ , 

/8 < 2 r pi(t)w{t)^dt < 2 r pi{t)w{t)t^^-'^+'dt 

So as in /s, 

< cie-^^" I pl{t)w{t)t^^-^^+Ut < cie-^^"^^ (30) 

That is /i - 2-^, and /a < cie '=2nv_!£z^ which proves the second hne of the 
lemma. 

For the proof of the Theorem we need the following 
Lemma 2 /// e C^, then by the notation of (13) 

\\YnA\^0{l) (31) 

Proof: 

At first let ^ <\x\< La„. Here, by W 5.5, - ^. Using that SJ£l 

2 — I I — " I J i a 1 Q (x) 71 ^ X 

is even, and it is increasing on when A > 2, we have 

Q {x) ^ w(xk) Xk 
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For the last equality cf. [12] (39). 

When |a;| > La„ > |a;fe|, then A>2 yields that g^^^f^ < 1- That is 

" 72 ^ ^ 

r„,^(x)<cu;(x)^^p- = 0(l), (33) 
fc=i ^^^''^ 

as in the previous case. 
Let d < < 

Yn,i,{x) = w{x) (O-Si + Es 

k k 
|xfc|<2|a;| |xfc|>2|a;| 

Because as in previously, when \xk\ < 2\x\, , ^^''^ < c(i3)' 

El - 0(1) (34) 

Y.^^w{x)vl{x)-^^ J2 ^(^,)^,(?l4^)2p;2(^^,) 

SFI^lk^T) ^ ('^f 4.11, 4.17) 

S2<cu.(x)p^(x)-^^ Y: -p^Ax, 
Q (X) n ^ Xk{x - XkY 



< cw{x)p^{x) 2n'f \ ~ dx<c — - 

x^Q (X) n X n xQ [x) 

Here we used that ■w{x)pn{x) < for |x| < cf. |L 4.6. 

Finally let \x\ < d. Let us remark, that is between two constants on 
[l,d]. 

Y^Mx)<cw{x) V iM + c^x) V iM^^ = E3 + E4 
As previously, 

Eg = 0(1) (36) 
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Similarly to the estimation of E2, 

E4<c^ E -p^Ax,<^^(^+Q{ar.))=0{l) (37) 

Proof (of the Theorem): Because polynomials are obviously in C^, ac- 
cording to the Banach-Steinhaus theorem, the previous lemma ensures the re- 
sult. 
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